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ABSTRACT 

The  stability  of  a  fluid  contained  in  a  circular  cylindrical  tank  with  a  flat,  flexible 
bottom  under  a  periodic  axial  excitation  is  studied.  A  variational  approach  is 
formulated.  An  approximate  solution  results  in  a  pair  of  coupled  ordinary  differ¬ 
ential  equations  with  periodic  coefficients.  A  method  of  handling  the  stability  of 
the  solutions  of  such  a  system  of  equations  is  presented.  Numerical  results  are 
discussed. 
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English  Symbols 


Constants 
Membrane  number 
Bond  number 
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g(t) 

g0 

gl 

G(r) 

h 


H 


L2’  L3 

»  M_2  *  M3 

N 


A  n 


CJ 


P 
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(r  ,  e  ,  z) 
(R  ,  0  ,  Z) 


Pi 

Amplitude  of  the  n  ‘  sloshing  mode 

Potential  of  the  solid-liquid-gas  interface 

Potential  of  the  edge  load  acting  on  the  rim  of  the  tank 
bottom 

Gravitational  acceleration,  time  dependent 
Mean  local  gravitational  acceleration 
Amplitude  of  the  imposed  axial  acceleration 
Nondimensional  gravitational  acceleration 
Membrane  thickness 
Nondimensional  free  surface  shape 
Functionals 

Bessel  functions  of  first  kind 
n^1  root  of  the  equation  (kQ)  =  0 
Depth  of  liquid 

Nondimensional  depth  of  liquid 
Pressure  energy  in  nondimensional  form 
Lagrangians  in  nondimensional  form 
Matrix 

Midplane  stress  resultant 

n^  sloshing  frequency  for  rigid  tank 

Pressure 

Nondimensional  pressure 
Radius  of  the  tank 
Cylindrical  coordinates 
Nondimensional  cylindrical  coordinates 
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id 


.  V  mm  -mtaanttrn  — 


mm 


Surfaces 


w 


W 

z 

Z 

Greek  Symbols 
a 


a 

r 


0 

$ 

CO 

n 


Time 

Modal  matrix 

Transverse  deflection  of  membrane 
Nondimensional  transverse  deflection  of  membrane 
Vertical  coordinate 
Nondimensional  vertical  coordinate 

Nondimensional  axial  acceleration 
Constants 

Boundary  curves  of  respectively 

Constants 

Free  surface  shape 

Azimuthal  coordinate 

Mass  ratio 

Constants 

Constants 

Density  of  membrane 
Density  of  liquid 
Surface  tension 
Nondimensional  time 
Velocity  potential 
Nondimensional  velocity  potential 
Forcing  frequency 
Nondimensional  frequency 
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,  82  (n  =  m  ,  ,  N) 


V2  ,  V2 


Eigenvalues  of  where  M  ,  ,  N  indicate  the 

fluid  modes  which  are  chosen  for  the  approximate 
solution 

Frequency  parameter  for  the  membrane 
Frequency  parameter  for  surface  tension 
Two-dimensional  Laplace  operator 

/  i  „2  -=,2  i  a  a  ^  i  a2  .  ..  .  .  , 

\^V  =  V  =  r  aT  r3~  +  ^  ^2  in  cylindrical 
coordinates^ 


V,  V 


Three-dimensional  gradient  operator 


V  =  v  =  e*  — ■  +  —  —  i  e*  — 

r  9r  r  30  ez  3z 


in  cylindrical  coordinates;  for  the  free  surface  and 


the  membrane,  ignore  the  e z  —  termj 


4 


THE  EFFECT  OF  WALL  ELASTICITY  AND  SURFACE 
TENSION  ON  THE  FORCED  OSCILLATIONS  OF  A 
LIQUID  IN  A  CYLINDRICAL  CONTAINER 


% 


1  INTRODUCTION 

Dynamics  of  large  liquid-fuel  rockets  naturally  involve  the  motion  of  a  liquid  in  a 
flexible  container.  The  symmetric  modes  of  the  fluid  motion,  which  influences  the 
pressure  at  the  tank  bottom,  and  therefore  influences  the  pressure  in  the  pump  and 
in  the  combustion  chamber,  as  well  as  thrust  and  rocket  acceleration,  have  an  impor¬ 
tant  effect  on  the  structural  dynamics  of  a  rocket.  In  some  instances  the  longitudinal 
oscillations  were  so  serious  as  to  affect  the  safety  of  the  vehicle.  For  this  reason, 
the  analysis  of  the  forced  oscillations  of  the  liquid  container  is  important. 

At  ground  level,  perhaps  the  effects  on  fuel  sloshing  of  the  flexibility  of  the  tank  wall 
and  the  surface  tension  of  the  free  surface  are  negligible.  At  reduced  gravity  condi¬ 
tions,  these  effects  become  more  evident.  It  is  the  purpose  of  this  article  to  evaluate 
the  effects  of  tank  flexibility  and  surface  tension  on  the  stability  of  liquid  motion  in  the 
symmetric  modes. 

Sloshing  of  liquids  has  been  studied  by  many  authors.  Although  most  of  them  considered 
rigid  containers  (Ref.  1),  Miles  (Ref.  2)  considered  bending  modes  of  a  flexible  con¬ 
tainer,  and  Bleich  (Ref.  3)  investigated  the  longitudinal  modes.  Recently,  Bhuta  and 
Koval  (Refs.  4  and  5)  studied  the  coupled  oscillations  of  a  liquid  in  a  tank  with  a  flexible 
bottom.  They  defined  the  normal  modes  of  the  system,  and  treated  the  orthogonality 
and  expansion  theorems.  Bhuta  and  Yeh  (Refs.  6  and  7)  considered  the  problem  of 
arbitrarily  assigned  velocity  distribution  on  the  tank  bottom. 

On  the  other  hand,  there  is  substantial  literature  about  the  influence  of  surface  tension 
on  sloshing,  e.g. ,  Yeh's  bibliography  (Ref.  8)  and  papers  by  Bond  and  Newton  (Ref.  9) 
and  Reynolds  (Ref.  10).  Most  of  these  studies,  however,  are  concerned  with  free 
oscillations.  Very  little  has  been  done  about  the  influence  of  surface  tension  on  forced 
oscillations,  and  no  work  seems  to  have  been  done  on  coupling  with  the  flexibility  of 
the  tank. 
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In  the  present  paper,  a  circular  tank  with  a  flexible  bottom  under  vertical  periodic 
excitation  is  studied.  The  problem  is  first  formulated  in  the  form  of  differential 
equations  and  then  in  the  form  of  a  variational  principle.  An  approximate  solution  is 
presented,  which  results  in  a  pair  of  coupled  ordinary  differential  equations  with 
periodic  coefficients.  The  stability  of  the  solutions  of  these  equations  is  discussed. 

2  STATEMENT  OF  THE  PROBLEM 


A  circular  cylindrical  container  with  rigid  side  walls  and  a  flat,  flexible,  bottom 
contains  a  liquid  with  a  free  surface.  The  tank  walls  are  subjected  to  an  oscillatory 
axial  acceleration,  in  addition  to  a  constant  mean-local-gravitational  acceleration 
directed  along  the  axis  of  the  cylinder.  Above  the  liquid  surface  is  a  gas  with  con¬ 
stant  pressure.  No  external  force  acts  underneath  the  tank  bottom.  The  situation  is 
pictured  in  Fig.  1.  The  problem  is  to  determine  the  motion  of  the  liquid  and,  in 
particular,  its  stability. 

The  fluid  properties,  including  the  surface  tension,  are  assumed  to  be  uniform, 
constant,  incompressible,  and  inviscid. 


The  mean  free  surface  of  the  liquid  is  assumed  to  be  a  plane  perpendicular  to  the 
cylinder  axis.  In  low-gravity  and  finite  surface  tension,  one  may  have  to  consider 
a  curved  mean  free  surface.  The  governing  criterion  is  the  Bond  number  defined 
below.  In  this  paper,  we  assume  that  the 
Bond  number  is  sufficiently  large  so  that 
the  free  surface  is  approximately  a  plane. 

The  case  of  low  Bond  number  is  discussed 
later. 


MIAN  ACCELERATION 


OSCItlATORY 
ACCfURATION 
o  COS 


MEAN  TREE  SURTACE  S 


As  a  further  simplification,  we  assume  that 
the  deviation  from  the  static  equilibrium  con¬ 
dition  is  small,  so  that  the  deflections  of  the 
free  surface  and  of  the  tank  bottom ,  the  fluid 
velocity,  and  hence  the  velocity  potential, 


RIGID  — 
WALL  S? 


FLEXIBLE  BOTTOMS. 


Fig.  1  Geometry  of  the  Problem  in 
Nondimensionol  Variables 
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may  be  considered  infinitesimal  quantities  of  the  first  order.  Under  this  assumption, 
all  the  equations  can  be  linearized,  and  the  mathematical  problem  is  relatively  simple. 
A  number  of  interesting  nonlinear  problems  are  ruled  out  by  this  assumption.  But, 
as  an  investigation  of  the  initial  tendency  toward  instability,  the  linearized  theory 
should  be  adequate. 

3  MATHEMATICAL  FORMULATION 

Consider  a  quantity  of  inviscid  liquid  situated  in  a  cylindrical  container  of  radius  r^ 
as  is  shown  in  Fig.  1.  The  cylindrical  polar  coordinate  system  is  chosen  so  that  the 
+  z  direction  is  directed  upward  away  from  the  liquid,  the  zero  on  this  axis  being  fixed 
on  the  mean  free  surface.  If  the  fluid  is  assumed  inviscid  and  incompressible,  and  the 
motion  irrotational,  the  equation  of  continuity  may  be  expressed  in  terms  of  the  velocity 
potential  <p , 


(i 


13  a  ,i  a2  ,  a2  i  ,  =  n 
3r  t2»f  az2'* 


and  the  velocity  components  u  ,  v  ,  w  are 


u  =  <t>r  ’  v  =  r  *0  ’  w  = 


The  usual  subscript  notation  is  used  to  denote  partial  differentiation. 


The  kinematic  conditions  at  the  tank  walls  and  the  free  surface  are 


u  =  fj  =  0  on  r  =  ro 


—  90  3w  „ 

W  =  8?  =  8t  °n  2  = 


w  —  — ^  —  ^2  on  z  —  0 
9z  at 


(1) 

(2) 

(3) 

(4) 

(5) 
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where  w  denotes  the  deflection  of  the  tank  bottom ,  and  t)  denotes  the  deflection  of 
the  free  surface,  both  positive  in  the  +  z  direction,  and  both  assumed  to  be 
infinitesimal. 

Since  the  motion  is  irrotational,  Bernoulli's  equation  is  satisfied  throughout  the  liquid 
domain.  In  particular,  at  the  free  surface,  we  have 

-  f(V<*>)2  -  g(t)r,  -  <f>t  +  c(t)  (6a) 

g  (t)  =  g0  +  gjCos  wt  (6b) 

where 

c  (t)  =  arbitrary  function  of  time 

gQ  =  mean  local  gravitational  acceleration 

g^cos  u)t  =  imposed  axial  acceleration* 
p  =  pressure  just  inside  the  interface 

The  pressure  p  is  related  to  the  pressure  just  outside  the  liquid,  p^  ,  by  the 
relation 


Taken  as  positive  if  they 
are  directed  toward  the 
tank  bottom  (along  the 
-z  direction) 


pG 


-  p  =  oK 


(7) 


where 

a  =  surface  tension 
K  =  total  curvature  of  the  free  surfaces 


♦Here  we  just  write  out  a  special  form  of  imposed  axial  acceleration.  The  method 
developed  later  can  be  applied  to  a  general  periodic  imposed  axial  acceleration. 
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In  linearized  form,  under  the  assumptions  that  77/r Q  «  1  and  |grad  tj|  <  1 ,  we 
have 


K 


_1_  8JJ 

r28* 


(8) 


If  the  pressure  of  the  gas  pp  is  a  constant,  then  without  loss  of  generality  we  may  set 
PG  =  0 .  The  function  c  (t)  can  be  absorbed  in  .  We  can  also  neglect  |  V0|  in 
Eq.  (6)  and  evaluate  on  the  surface  z  =  0  under  the  scheme  of  linearization.  Thus 
we  obtain  the  linearized  free  surface  condition, 

+  g(t )T)  (9) 

p0  '  l'z  =  0 


Similarly,  Bernoulli's  equation  gives  the  pressure  on  top  of  the  bottom  wall 


(P)z  =  -£+w  ^0 


-  (tf>t)  -  g(t)  (-f  +  W) 

z  a 


(10) 


No  other  forces  are  assumed  to  be  acting  on  the  tank  bottom .  If  the  tank  bottom  is 
very  thin  and  is  prestressed  so  that  it  behaves  like  a  membrane,  then  the  equation  of 
motion  of  the  bottom  is 


NfV2w  =  ph 


a2w 


at 


f  +  *<*> 


+  (p)  . 

VF/z  -  -J0  + W 


(11) 


where 

Nr  =  tensile  stress  resultant  in  the  tank  bottom  (assumed  to  be  a  constant) 
p  =  density  of  the  tank  bottom  material 
h  =  tank  bottom  wall  thickness 
ph  -  mass  per  unit  area  of  the  tank  bottom 
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A  combination  of  Eqs.  (10)  and  (11)  gives  the  linearized  equation  of  motion  of  the 
elastic  bottom  as  a  membrane: 

2 

NrV2w  =  ph^  -  P0 g ( t ) w  -  pQ|||j  ^  +  (ph  +  P0*)  g(t)  (12) 

In  reality,  a  tank  with  flat  bottom  develops  both  bending  and  stretching  stresses  under 

fluid  pressure.  Equation  (12)  is  a  good  approximation  only  if  a  membrane  tension  is 
built  in  at  the  edges  by  stretching  the  bottom  onto  a  rigid  cyi  ader  before  the  two  are 
welded  together. 

It  is  necessary  to  specify  the  boundary  conditions  for  77  and  w  at  the  edge  r  =  r^  . 
We  choose 

w  =  0  when  r  -  rQ  (13) 

|^  =  0  when  r  =  rQ  (14) 

The  last  condition  is  a  special  case  of  zero  capillary-hysteresis.  It  is  consistent  with 
the  simplifying  assumption  that  the  undisturbed  free  surface  is  a  plane  z  =  0  .  In  a 
very-low-gravity  condition,  the  mean  free  surface  is  curved,  and  Eq.  (14)  should  be 
replaced  by  the  condition  8  77/8  r =  yrj  at  the  wall  where  y  is  a  physical  constant. 

These  equations  define  the  linear,  inviscid  problem  of  sloshing  under  appropriate 
initial  or  periodicity  conditions. 

4  DIMENSIONLESS  EQUATIONS 

Taking  the  radius  of  the  cylinder  rQ  as  the  characteristic  length,  the  gravitational 
acceleration  gQ  as  the  characteristic  acceleration,  and  u>  as  the  characteristic 
frequency,  we  define  the  dimensionless  variables  as  follows: 
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R 


_r 

r 


0 


Z  =  -5. 
r0 


*r 


L 


$ 


_£ 

r 


0 


= 

2 

"ro 


W  =  — 


0 


G(r)  = 


g0  +  cos  w t 
g0 


T  —  Cjt 


H  = 


a  = 


P  = 


JL 

ro 


h 

h 


2  2  i 

%  "  r0  / 


(15) 


We  define  the  dimensionless  parameters  as  follows: 

2 

Bond  number  =  B^  =  pQ  gQ  rQ  /a 


Membrane  number  =  B„.  =  p„g.r“/N 

M  0  0  0  r 


2  3  2 

Frequency  parameter  for  surface  tension  =  =  pQ  r^u?  /a 


(16) 


2  3  2 

Frequency  parameter  for  the  membrane  =  =  pQ  rQ  w  /Nf 


Mass  ratio  =  \  = 


po  ro 


and  the  operator 


V2  =  1/R  8/8R  (R  9/3R)  +  l/R2  92/802  =  r2  V2 


(17) 


Then  the  equations  become 


+  v2  )  *  =  o 

8Z2  / 


(18) 
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V2  H 

^ \  -  B  G  H  =  0 

*T/z=o  a 

(19) 

V2  W  -  A  fl2  82f+fl2  (»♦ 
M  ar2  M  \9T 

■u 

+  BmG(t)W  -  (X  +  L)BmG(t)  =  0 

(20) 

with  the  boundary  conditions 

S|  =  0 

aR 

on 

R  =  1 

(21) 

a*  _  aw  „„ 
az  3t 

Z  =  -L 

(22) 

a*  an  _ 
az  3t 

Z  =  0 

(23) 

W  =  0 

on 

R  =  1 

(24) 

and 

M  =  o 
aR  u 

on 

R  =  1  (assuming  y  =  0 ) 

(25) 

Equations  (18)  through  (25)  show  that  the  problem  of  sloshing  depends  on  the  parameters 

2  2 

£2  ,  B  ,  a  ,  £2*.  ,  B,.  ,  X  ,  and  L . 

<7  a  M  M 

These  dimensionless  parameters  are  not  all  independent;  since 

°M  "  N"  ““  <26> 
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therefore 


2  2 

However,  we  retain  the  sets  of  symbols  and  ,  BM  because  these  two 

pairs  of  parameters  are  not  both  likely  to  be  important.  The  conditions  are: 

2 

*  ,  B^  —  0  if  the  tank  bottom  is  rigid 
2 

•  0  ,  B  —  oo  if  the  surface  tension  has  no  effect 

cr  1  a 

5  DISCUSSION  OF  ANALYTICAL  SOLUTIONS 

Consider  symmetric  modes  of  motion  in  which  $  ,  H  ,  W  are  independent  of  the  angular 
coordinates  0  .  A  solution  of  Eq.  (18)  may  be  posed  as 


^  cosh  k  Z  sinh  k  Z 

do<T>z*VT>  +  2Jo(k'>R)  VT>5snrL*VI>;sanr-L 

“  n  n 

n  =  l 


Then  Eqs.  (22)  and  (23)  give 


=  dA<T)  + 


OO 

2“ 
L.  « 


lS,J.<"n11> 

cosh  k  L 
n 


OO 

vr> +  f(*> +  -  vt>]  vs.R> 
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Both  Eqs.  (21)  and  (25)  are  satisfied  if  the  k^'s  are  the  roots  of  the  equation 

JL(kn)  =  0  n  =  1  ,  2  .  3  , .  (31) 

Eq.  (24)  is  satisfied  by  taking 


OO 


:<T>+  2  kn[dn<T>  -cn<T> 


'0 


<k„> 


=  0 


n=l 


(32) 


and 


R)  =  (X  +  L) 


1 


J0  R  > 
Jo(\f%> 


Here,  f(H)  is  the  static  deflection  of  the  membrane.  We  assume  that,  if  BM  is 
positive ,  is  less  than  the  first  root  of  Jq  (x)  =  0  ,  namely,  2.4048.  To  satisfy 

Eqs.  (19)  and  (20),  we  substitute  4>  ,  H  ,  W  from  Eqs.  (28)  through  (30),  collect 
terms,  and  represent  the  lefthand  side  as  a  Fourier-Bessel  series  in  (k^  R).  Since 

the  series  vanishes,  every  coefficient  of  (k^  R )  ,  n  =  1 ,  2  , - ,  must  vanish. 

Thus,  from  Eq.  (19)  we  obtain  the  necessary  conditions 

nJc0(T)  +  BffG(r)d0(r)  =  0  (33a) 


and 


(33b) 


Multiplying  Eq.  (20)  by  R  and  integrating  with  respect  to  R  from  0  to  R  ,  we  get 
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aw 


R 

/ 


R 


Rat  +  bmg  I  RW(R>T>dR  -  *oM 


5/ 


W(R,r)  RdR 


R 


lA-i  L)  n  PR2 
2  BM GR 


+  4  j  §7 


Z=-L 


RdR  =  0  (34) 


2 

Any  function  W(R,t)  of  class  C  in  the  closed  interval  0  to  1  for  R  satisfying 
the  above  equation  will  satisfy  Eq.  (20).  A  substitution  of  Eqs.  (28),  (29),  and  (30)  into 
it  gives 


OO 

2 

n=l 


d0R 


OO 


2  + 

n=l 


-cn)Ji(k  R) 


-x  a 


M 


doR 


OO 


n=  1 


+  0 


\  o  +  c' 


0ni  V  Jl(knR) 


M 


R  + 


n=  1 


n 


( c  coth  k  L  -  d  tanh  k  L ) 
'  n  n  n  n  7 


a 


(X  +  L)^ 

vv%>  Jl 


(V%R)  = 


0 


(35) 


By  expanding  R  and  (\/bm  R )  in  terms  of  (  k^  R) ,  we  can  collect  the  coefficients 
of  (  k^  R )  and  set  them  equal  to  zero,  to  obtain 


IT  lSa<“*h  k„  L  *  *  V  '  4n«ta"h  ^  L  *  X  V I  *  (kn  -  bm  gK  ‘  V 
n 


M 


knVkn> 


bm  G 

( X  +  L )  dQ  +  Cq  +  — ^ —  d. 


2  "0 
M 


2tt(X+  L)VB^Jl(^)  kn 


kn  J0(  kn>  J0(VBM)  Bm  -  kR 


(36) 
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Now,  if  we  want  to  truncate  the  infinite  series  in  Eqs.  (28)  through  (30)  by  taking 

n  =  m, . ,  N  ,  we  see  that  Eqs.  (32),  (33),  and  (36)  always  involve  2(N  -  m  +  2) 

unknowns  and  2(N  -  m  +  2)  equations,  which  in  general  have  solutions. 


Of  course,  we  can  do  the  same  thing  for  Eq.  (19)  to  obtain 

o  .  /  9  v  sinh  k  L 

ic n  +  K  +  Ba  G)  Vtanh  kn  L>dn  =  XTFT  ( 


W  rvo 


n2  C.  +  G  dQ)  (37) 


However,  since  we  have  assumed  8H/8R  =  0  at  R  =  1  ,  the  series  of  Eq.  (29),  after 
twice  term-by-term  differentiation,  is  still  convergent.  Thus  Eq.  (37)  coincides  with 
Eq.  (33). 

We  shall  now  discuss  the  solution  for  Eqs.  (32),  (33),  and  (36).  For  the  case  of  free 

vibration,  i.  e.  G  =  1  or  a  =  0 ,  solutions  for  c  ,  d  can  be  obtained  rather  easily. 

n  n 

If  we  let 


c 


n 


ittT 

e 


and 


(38) 


d  =  D  e 
n  n 


i  Qt 


then,  from  Eq.  (33),  we  obtain 


C 


0 


B 


2  2 
ST  U 

<7 


n 

~0 


and 


(39) 


+  B 
_ a_ 

n2  n 

a 


knL 


n 


\ 
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_  .  ~  tea*  >  Miami 


and  from  Eqs.  (36)  and  (39), 


2 

J0(kn)D0 


P  D 
rn  n 


where 


R2  P2  (coth  k  L  +  X  k  ) 
M  *n  '  n  n ' 


R2  fi2(tanh  k  L  +  X  k  ) 
M  '  n  n' 


)(5 


Substituting  into  Eq.  (30),  we  have 


eiOTD0  1  + 


00 


n2  n2(L +  x) 


p^\Jrt(k  R) 

?  w  +  ,<RI 


Now  we  can  easily  see  that  the  coefficient  of  J  (  k  R )  in  the  above  series  is  of  order 

o '  n 

l/(k3/2)  for  large  n  ;  so  the  series  is  actually  divergent  after  twice  term-by-term 

2 

differentiation  with  respect  to  R.  At  R  =  1  ,  W  =  0  ,  we  have  D  ^  0  for  fit,  &  0 

o  M 

2 

(S^M  s  0  corresponds  to  rigid  tank);  therefore, 


f(0  )  =  1  + 


WL*"\,2  2 V  „ 

- T - (°  ‘  pn)  ‘  0 
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2 

This  equation  will  determine  the  eigenvalue  ft  .  It  can  be  shown  that,  for  Pn  -  ^  t 

n 

^  ^0  ,  all  the  roots  ft  are  real,  and  no  double  roots  exist.  In  the  case  — 0,  some 
ft!s  can  be  obtained  asymptotically  in  a  rather  simple  fashion.  Let 

°2  *  -  A  (>  -  4  -  S2»M  -  «J<4  ...  )  («) 


Then  we  have 

p£  ~  [Pj j  ^  )  “  P^  ~  ^  “  ^2  *  *  *  )  ^  + 

**l8l0Mt(42  +  4l)®»l*  ■•I 


=  [dj  X|  -  Pj  (coth  kj  L  -  tanh  kfL)  J  £2^ 

+  [  (61  +  ^2)  *«  "  p«  61 (tanh  k«  L  +  Xkf  >)  nM 

+  [  («J  +  26x  62  +  63)  Xf  -  Pf  62  (tanh  kf  L  +  X  k,)]  £2^  +  0  (l^)  (45) 

where 

X|  *  k.(k?  -  bm)  «46' 

Equation  (43)  becomes 


1  - 


Pf  ( L  +  X)(61  +  62  £2^  +  ...  )ff 


2 

'M 


6lxt  "  P{  (°°th  kf  L  -  48,111  kf  L>  +  [  (6i  +  62 )  x«  “  Pjf  (tanh  kf  L  +  X  1^)]  £2^  +  . . . 


nMp«(L  +  X) 


00 

I 


n  =  l 


kl  "  bm 

n  M 


k?‘BM 


+  0 


(<4)  ■ 


(47) 
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is.  istmt* 


In  order  that  the  above  equation  be  valid  as  —  0  , 


we  must  have 


6 


1 


coth  L  -  tanh  k|L 
Xjt  ” 


2 

pi 


and,  since 


(48) 


(49) 


Eq.  (43)  becomes 

_ 

(S1  +  s2 )xl  '  P|' 51<tanh  kl  L  +  +  |xj(61  +  261*2  +  6a)  ‘  P|2l52<tanhk«L  *  Xk«>)  nM  +  •" 


-  «mpi  <L  +  M 


0 


(50) 


Expanding  the  lefthand  side  in  power  series  of  flM ,  and  putting  the  coefficient  of 
2n 

equal  to  zero,  we  get 


Pt  (L  +  X) 


(61  +  ^2)  *f  ‘  Pf6l<tanh  k|L+  Xkf> 


=  0 


(51) 


and 


Pj  (L  +  X)  62  -  (62  +  26^  +  63)  x{  +  p^<52(tanh  kf  L  +  Xkf) 


[p2(l  +  X)| 


(V=S) 


B 


M 


2V^Ji (V®m)  kf  ”  bm j 


6=0 
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f>  Min 


Then 


*2  = 


v  ]vi  - 

*t 


and 


(52) 


6„  = 


_1_ 

Xf 


*  2 
62PiVl 


* 

-(^*2V2)x,*|'’?<L+H 


2>/®m  ji  (yfhi ) 


kf-»M 


Therefore,  after  some  rearrangement,  we  get 


2 

pf 


i  - 


h  _  iy*ipi' 


A 


f 


-  4Vi 

k  -i  A I 


where 


‘I 


(53) 


a  ii 

f  "  fi2 
*  M 


(54) 


From  Eq.  (53)  we  make  two  interesting  observations: 

•  The  term  is  the  f  ^  nondimensional  natural  sloshing  frequency  of  the 
liquid  in  a  rigid  tank.  Since  (p^/^  )  >  1 ,  we  conclude  that  <  Pf  . 

•  The  elastic  effect  is  at  least  of  order 


a 


M 


sinh  2kJf  L 


=  2 


T.  ("i  -  "«)] 
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Because  of  the  presence  of  the  elastic  bottom,  the  sloshing  mode  shapes  ire  no  longer 
as  simple  as  those  in  a  rigid  tank;  but,  in  rigid  tanks,  the  influence  of  cross  coupling 
of  different  fluid  modes  on  the  natural  frequency  is  of  order  . 

o 

The  other  limiting  case  is:  Pq  ,  the  density  of  the  fluid,  tends  to  zero  while  B^,  , 
are  bounded  away  from  zero;  then  Eq„  (43)  will  tend  to  the  free-vibration  frequency 
equation  for  a  circular  membrane  -  i. e. , 


£ 


°2>  *  >  +  I 


kV 


'  k2  Q2  -  k2 

n=  1  Kn  Kn 


kfi  J0  (kfi) 

^J1(k  fi) 


(55) 


or 


Jo(kQ)  =  0 


where 


2 


6  STABILITY  OF  THE  SOLUTION 


To  study  the  stability  of  the  solution  (Refs.  11  and  12),  we  shall  consider  the  following 
more  general  system  of  equations:* 


••  .  A  i2nt  ,  v  a  i2nt 

I  A2ne  A2n-le  ^2  =  0 


n=  -  oo 


n=  -  co 


(56) 


♦This  is  an  extension  of  Hill's  method  (see  Ref.  12,  p.  4.13)  to  a  system  of  two  equations. 


where  j3q  *  Aq  ,  Aq  >  0  ,  >  0  ,  A  ^  =  B,  =  0  ,  ^  0n  (where  /?n  =  | j  A  |  cr 

n=-<» 

|Bn| )  is  an  absolutely  convergent  series.  Equations  (56)  are  invariant  when  t  is 
changed  to  t  +  it ;  therefore  if  y  ( t )  is  a  solution  of  Eqs.  (58),  y(t+7r)  is  also  a  solu¬ 
tion.  By  the  Floquet  theorem  (Ref.  13),  Eqs.  (56)  have  solutions  of  the  following  form: 

y  (t)  =  e?t#>(t) 


where 


1  = 


*2 


The  term  ip(t)  is  periodic  function  mod(?r) .  If  Re  {  >  0  ,  y as  t-*  an  unbounded 
solution  exists,  which  is  said  to  be  unstable.  For  a  periodic  solution  mod  (7r)  to  exist, 
lm£  must  be  equal  to  an  integer,  whereas  Re  4  =  0 

Let  us  assume  a  solution  of  the  following  form: 


...  \  i2nt 

yl(t)=e  Z  r2ne 


n=  -  co 


...  $t  V  i2nt 

y2<l)  =  e  >  y2n  + 1  e 


n=  -  00 


2  ^  2 


where  ^  n  «  HI  n*-  ^2n+l  are  a^so^u^e^y  convergent  series. 


n-  -  00  n~  -  °o 
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Substituting  into  Eqs.  (56),  we  get 


00 


oo 


OO 


2  Vau  +  e, 2  ^  A2me2»‘*  2 


Y2ne 


(2ni  +  £)t 


n=  -  « 


m=  - » 


n=  -  oo 


OO  OO 

5  2 


*2n+l  e 


n=  -» 


(2ni+4)t  =  Q 
(57) 

2  W2"i+<>2e|2"i*S|t*  ^  B2m  e2mi‘  2  ^2n+ 1  ®,2n‘ 


m=  - 


OO 


£)t 


n=  -  oo 


m=  -  oo 


n=  -  oo 


2  B2m+1  e2mit  2  ^‘(2nln,‘  '  • 


m=  -  °o 


n=  -oo 


On  rearranging  the  terms  of  the  absolute  convergent  series,  and  equating  the  coefficients 
of  +  ^  to  zero,  we  obtain 


oo 


_Y  l-il.--.2nr  +  y  -m 

2nAn-4n2  Z  A.  -  4n2  ^2n_m 
0  m=  -  oo  0 


-  0 


oo 


(i£  -  2n )  ,  \ 


B 


m 


'V2n+1TT^+  Z  B0-4n2  2n+1'm 
0  m=  -  °o 


=  0 


(58) 


2  o 

provided  that  AQ  -  4n  *  0  ,  B0  -  4n^  *  o .  The  divisors  AQ  -  4n2  and  BQ  -  4n2  are 
introduced  in  order  to  make  an  infinite  determinant,  which  will  be  formed  below,  to  be 
convergent. 
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Equations  (58)  are  a  set  of  homogeneous  equations.  For  to  have  nontrivial  solutions, 
the  determinant  formed  by  the  coefficients  of  the  equations  must  vanish.  Call  this 
determinant  A  ( 4  ) ;  then 


A(4 )  =  \ot„\  =  0 


(59) 


where 


a 


2m,  2m 


Aq  -  (4  -  2m ) 
Ao  -  4“2 


a 


2m+l  ,  2m +  1 


B0  -  (i?  -  2m)i 
B0  -  4m2 


tt2m,n 


2m -n 

Ao  -  4“ 


for  2ra  -  n  *  0 


a 


2m+l ,  n 


B 


2m  + 1  -n 


B0  -  4m 


2 


for  2m  +  1  -  n  =  *  0 


m ,  n 


=  0  ,  ±1  ,  ±2,.. 


We  consider  another  infinite  determinant  (4)  =  |/3y|  where 


B  =  1 
Mm ,  m 


(60a) 


a 


2m  ,  n 


2m  -n 


2m,  n  0  »  ...  0  ,2 

2m, 2m  AQ  -  (4  -  2m) 


for  2m  -  n  *  0 


(60b) 


a 


2m  + 1 1  n 


B 


2m+ 1  -  n 


2m+ 1 ,  n  a2m+ 1 ,2m+ 1  BQ  -  (i£  -  2m)1 


for  2m+l-n^0  (60c) 
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Since 


77  /*  =  1 

I  J  rm,m 
m=  - 00 


m,n=  - 00 
m  *  n 


1  r  1 


converges,  provided  4  does  not  have  such  a  value  that  one  of  the  denominators  of 
Pmn(*£)  vanishes.  Thus,  the  infinite  determinant  A^  ( i£  )  is  absolutely  convergent. 
Then  (Ref.  14), 


A(4)  =  Ax(i« 


lim 
m— 00 


[a0  -  (i€  -  2n)2 

[bq  -  (4  -  2n)2 

(a0  -  4n2) 

(Bo  -  4”2) 

We  note  some  interesting  properties  of  A^  ( i£  ) :  (1)  A^  ( 4  )  is  a  meromorphic  function 
of  4  and  tends  to  1  as  Re  4  —  ±<*>  ;  (2)  A^  (4  )  is  a  periodic  function  of  4  with  period 
2i  .  If  we  form  another  function, 


F(4  )  =  \  (i| )  -  [cot|(i£  +  v^)  -  cot|(i$  - 

-  K2  |cot|(U  +  v7^)  +  cot f  (i?  -  v^)j 

-  K3|cot|{i4  +  -JW^)  -  cot | ( 4  -  /B^)j 

-  K4(cotf(4  +  vEJP  +  cot|(i4  - 


(62) 
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where  K/s  are  so  chosen  that  F(£)  has  no  poles  at  i£  =  ,  ±vT3^,  then,  since 

A^  ( 4  )  is  a  periodic  function  of  5  ,  it  follows  that  F  ( £  )  *  has  no  poles  at 

4  =  2n  ±  yfK^  ,  2n  ±  ,n=  ±  1  ,  ±  2  , - 

Thus,  F  (^  )  is  a  meromorphic  function  with  no  pole  on  the  entire  plane.  F  (£  )  is 
certainly  bounded,  therefore,  by  Liouville's  theorem,  F  ( M  must  be  a  constant,  say 
C .  As  ReJ  -  ±oc  ,  ^(i^)  =  1.  Therefore, 

C  =  1  +  2  ( K2  +  K4 )  i  as  Re  £  —  °o 

=  1  -  2  ( Kg  +  )  i  as  Re  £  —  -<» 

Hence,  Kg  +  Kg  =  0  ,  and  F  (£  )  =  1  for  all  £  .  Using  this  result  and  Eqs.  (61)  and 

(62),  we  get 


A(4)  = 


sin4  -  262  sin2  ^  +  63  sin  ^ 

sin2  (f  *£) sin2  (f  ^o) 


(63) 


where  6i  are  some  constants  relating  to  K^'s  ,  ,  and  BQ  .  Put  4  =  0,  1/2  and 

1  in  Eq.  (63)  and  we  get 


“i  ■  M0)Si»2(f^)si„2(f/ir0) 

262  -  1  *  sin2(|^)  sin2 1^(0)  -  a(l)|  (64) 


and 

263  =  |  +  sin2(fv^)  sin2(f  jwj)  |2a(|)  -  A(0)  -  A(l)j 
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f 


4 


In  a  special  case,  if  the  coefficients  of  Eq.  (56)  are  even  functions  of  t  ,  then  Eq.  (56) 
is  unchanged  when  we  change  t  to  -t  ;  we  see  that,  if  $  is  a  solution,  then  -£  is  also 
a  solution.  Therefore,  if  A(i£Q)  =  0  and  sin(i£07r)*  0  ,  then  6^  =  0  .  Therefore, 
when  we  want  to  find  the  roots  for  A(i£)  =  0  ,  we  always  have  6Q  sin  ( n.  4  )  =  0  and 
the  roots  of  Eq.  (63)  can  be  written  out  in  a  simple  form, 

sin2  ^  =  <52  *  "  61  (65) 

For  a  boimded  solution,  i.  e.  for  Re£  =  0,  we  must  have 

1  a  S2  *  "  61  £  0  (66) 

In  Eq.  (63),  by  putting  A(i£  )  =  0  ,  we  can  compute  £  »  and  determine  whether  this 
is  an  unbounded  solution  or  not.  Then,  from  Eq.  (58),  we  can  compute  y '  ,  and  obtain 
the  complete  solution  of  Eqs.  (56).  For  a  periodic  solution,  we  must  have  A(0)  =  0  or 
A  (1 )  =  0. 

If  a  periodic  solution  of  an  inhomogeneous  counterpart  of  Eqs.  (56)  is  considered,  e.  g. 

1  +  A(t)x  =  B(t) 

where  B  (t)  is  a  column  matrix  with  its  elements  as  periodic  function,  we  can  use 
Eq.  (57)  (with  non-zero  righthand  side)  by  putting  g  =  0  ,  whereas  Eqs.  (58)  become 
inhomogeneous.  If  A(0)  is  not  equal  to  0  ,  we  can  solve  for  y  uniquely.  If  A(0)  =  0 
we  are  on  the  boundary  where  Eqs.  (56)  have  an  unbounded  solution.  Therefore,  for 
such  an  inhomogeneous  equation  as  in  our  problem,  the  zone  of  instability  is  determined 
by  the  homogeneous  solution. 
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7  APPROXIMATE  SOLUTIONS 


We  choose 


*  =  d0(T>z  +  co<r)  + 


J  r  cosh  Z  .  sinhknZ 

^  Jo^knR^  I  cn ^  sinh  k_  L  dn^  cosh  k_  L 

n  =  m  L 


n 


n 


(67) 


N 


H  = 


VoV 


d  (t)  +  \  d  (r) - ,  .  T 

°  Z*  n  cosh  k- L 

n  =  m 


(68) 


n 


and 


N 

W  =  d0(T)  +  2  kn  ldn(T)  '  Cn(T))  Jo(knR>  +  f(R)  (69) 

n  =  m 


where 


d0<T> 


N 


'  5  kn  [  cn<T»  *  dn<T>  ] 


j  (lC  ) 
o'  n7 


n  =  m 


(70) 


These  equations  satisfy  Eqs.  (18),  (19),  and  (21)  through  (25).  Then  the  equations 

fia’6o(T)  +  BaG(T)do(T)  =  0  <71a) 

fl^cn(T)  +  kntanhknL  [  k^  +  G(t)  j  dn(T)  =  0  (71b) 

hold  for  n  =  m,m+l,...N.  Now,  N  -  (m  -  1 )  of  the  cn  ,  dn  are  left  arbitrary. 
To  determine  these  arbitrary  functions,  we  use  the  Ritz  method  by  substituting  4> , 

H  ,  and  W  of  Eqs.  (67)  through  (70)  into  the  variational  equation 


6 


(72) 


f  I  (*  ,  H,  W  ;  T)dr  =  0 
1 

The  Euler  equations  are  exactly  Eq.  (36)  for  N  =  m  . . .  N ,  with  cQ  ,  dQ  satisfying 
Eq.  (71a). 


For  free  vibration,  by  taking  a  single  term  m  =  N,  we  get  the  none- term  approximation" 


1  + 


=  0 


(73a) 


or 


fll 

ft2 


1  2 
2  Pm 


Am  Mm  - 

2  v  + 

v  p  m 
m  m 


//A 

fi  \ 

2  4  A 

y  m 

m  > 

111 

(..  2 

V 

..  .2 

^  Prv, 
mm 


m 


m  rm 


(73b) 


which  is  the  same  as  neglecting  all  the  terms  in  the  summation  of  Eq.  (43)  except  the  m**1 
term.  If  we  take  two  terms,  i.e. ,  N  =  m  +  1  ("two-term  approximation"),  we  get  the 
following  characteristic  equation  for  ft  : 


n~m 


(74) 


8  NUMERICAL  RESULTS 


Some  numerical  results  are  obtained  using  the  method  discussed  above.  In  Figs,  2 
through  5,  the  natural  frequencies  normalized  by  the  sloshing  frequency  p^  ,  which 
equals 


tanh 
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of  the  corresponding  rigid  tank,  versus  the  membrane  number  B.. ,  which  equals 

(c*  >  M. 

p  g  r*j/N  ,  were  plotted  for  parameters  B  ,  L,  A,  Since  ft  is  normalized 
o  o  o/  r  2  2  ®  2  1  \ 

by  ,  the  parameters  ft^  and  ftM  ,  which  equals  )  can  be  eliminated. 

The  computation  was  based  on  Eq.  (73)  (one-term  approximation),  with  m  =  1  ,  and 

on  Eq.  (43).  In  most  of  the  calculations,  the  one-term  approximation  and  exact  solution 

gave  almost  identical  results  for  the  lowest  frequency,  which  cannot  be  distinguished  in 

the  figures  shown. 


From  these  figures,  we  see  the  general  trend  clearly  that  frequencies  decrease  as 
Bm  increases.  In  Fig.  2,  the  Bond  number  is  large  (B^  =  10),  and  B^  does  not 
have  much  effect  on  the  lowest  frequency.  As  B^  decreases  (see  Figs.  3  and  4, 

B^  =  1 ,  0.1),  the  great  effect  of  B^  can  be  seen.  At  certain  ranges  of  B^  ,  the 
lowest  frequency  decreases  sharply  as  B^  increases;  and  the  range  is  approximately 
determined  by  the  ratio  of  B^.  to  B^  ,  and  the  depth  of  liquid  L.  This  relationship 
can  be  seen  clearly  from  Eq.  (43),  if  the  terms  in  the  summation  sig.  are  normalized 
in  a  slightly  different  form ,  i.e.,  for  B^  ,  B^  small, 


1  + 


kn(L  +  \) 

1  /«y  (k? +  %)  ki  tanhRi  l 

lpl/  (k2  +  bJ  k  tanh  k  L 

t  x  1  \  n  o!  n  n  J 

k2  '  bm 
n  M 


/ k ^  +  B  )  k-  tanh  k.  L 

_ ,  coth  k  L  +  Ale - «  Vo - r - +  r2  .  r  I  2 

n_1  pf  k  +  B  )  k  tanh  k  L  (k  +  tanh  k  L  UM  \p 

*1  \  n  (J/  n  n  \  n  cr /  n  '*n 


k„  (L  +  X) 


j  _,iL\2  kitanhkiL 

ku  tanh  k  L 
n  n 


n  =1  c°th  kfiL  + 


k1  tanh  k2  L 


B 


k l  tanh.  kn  L  +  tanh  kn  L  BM 
-L  n  n 


^\2  k^  tanh  k^  L 

P1  /  k^  tanh  k  L 
n  n 


-  1 


=  0 
(75) 
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Fig.  2  Frequencies  Computed  From  Eqs.  (73)  and  (43), 
With  Bond  Number  10 


Fig.  3  Frequencies  Computed  From  Eqs.  (73)  and  (43), 
With  Bond  Number  1.0 
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Apparently  the  roots  for  ft/p^  depend  on  B^/B^  only  for  fixed  L,  X.  Physically  it 
says  that,  at  low  Bond  number,  the  surface  tension  plays  an  important  role;  its  effect 
depends  mainly  on  B^/B^  (the  ratio  of  the  surface  tension  to  midplane  stress  resultant 
of  the  membrane).  The  argument  holds  for  any  N-term  approximation. 

In  Figs.  4  and  5  we  can  see  the  effect  of  the  depth  of  the  liquid  to  radius  ratio  L .  It 
causes  the  sharp  decrease  of  frequency  occurring  at  smaller  B^  for  larger  L . 

In  Fig.  6  (a  and  b),  the  stability  boundary  with  o:  (=  g1/gQ )  is  plotted  versus  forcing 
frequency  normalized  by  the  rigid  tank  first  sloshing  frequency,  according  to  Eq.  (A-9) 
in  the  Appendix,  with  i  =  j  =  1  for  one -term  approximation  only. 

9  CONCLUSIONS 

From  the  exact  solution  for  free  vibration,  we  make  four  conclusions  for  the  eigenvalue 

2 

ft  :  (1)  ft  is  always  real,  for  p^  ,  xn  >  0  ;  (2)  the  elastic  effect  lowers  the  natural 

vibration  frequency  ;  (3)  the  elastic  effect  on  each  sloshing  frequency  is  of  order 
2 

ft^/sinh  2k^  L  ;  (4)  because  of  the  presence  of  the  elastic  bottom,  coupling  of  different 
sloshing  modes  occurs,  and  the  effect  of  one  sloshing  mode  on  the  natural  frequency  of 
another  due  to  this  elastic  coupling  is  of  order  ft^  . 

Numerical  results  indicate  a  great  effect  of  surface  tension  on  the  natural  frequency  and 
stability  boundary  at  low  Bond  numbers.  If  B^  ,  B^  are  both  small  (  <1  ) ,  then  the 
frequency  depends  on  the  combination  of  parameters  B^/B^ . 
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Appendix  A 

APPROXIMATE  BOUNDARY  FOR  A  SYSTEM 
OF  DIFFERENTIAL  EQUATIONS  WITH  PERIODIC  COEFFICIENTS 
TO  HAVE  UNBOUNDED  SOLUTION 


Following  Struble's  (Ref.  15)  approximate  stability  boundary  for  Mathieu's  equation, 
we  extend  it  to  a  system  of  N  equations: 
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^  +  Dy  =  a  cos  t  Ay 


(A-l) 


where 


A  - 


(A-2) 


and  ft.'s  are  different  from  zero,  a  is  a  small  parameter, 
assume  that  y.’s  have  the  expansions 


a.,  are  constants.  We 
i] 


y.  =  E.(t)  cos(i2.t  -  6.  t)  +  a  y|1}  +  a2  y|2)  + .  (A-3) 

•  2  • 2  ••  •  • 

where  E.(t)  „  0.(t)  are  slowly  varying  parameters,  and  that  E.  ,  0.  ,  E.  ,  0.  E.  , 
0^  can  be  neglected  as  compared  with  E.  ,  0.  .  Substituting  (A-3)  into  (A-l),  we  get 


2J2.  [e.  el  +  0 (62 ,  E2. . . .) ]  cos  (fi.  t  -  6 i )  -  2£2.  (E.  +  . . . . )  sin  (fi.  t  -  9 t)  +  a  [y}1 )  +  0?  yj1  ■ )] 
n 

f-  ^  a..  E.  (cos [(1  +  n j)  t  '  0 j  J  +  cos ((1  -  n.)  t  +  0  J  |  +  0 (a2)  i  -  1,2, . ,n 

j=l 

(A-4) 
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For  ft.  *  1  ±  ft^  ,  we  may  solve  for  immediately  by  taking  E.  ,0.  to  be  constants. 

But  there  is  a  possibility  that  the  expansion  for  y^  breaks  down  because  of  zero  or 

small  divisors  for  the  solution  of  y^ .  *  In  order  to  avoid  this  difficulty,  we  can  use 

,0.  to  cancel  those  terms  that  cause  trouble.  We  shall  consider  the  case  that 

ft.  +  ft.  =  1  +  0(0?)  for  only  one  pair  of  specific  i ,  j . 
i  J 

In  this  case,  we  remove  the  troublesome  terms  in  (A -4)  by  putting 

2ft.  0.  E.  cos  (ft.  t  -  0. )  -  2ft.  E.  sin  (ft.  t  -  0. ) 

lii  'i  i '  ii  'i  i 

=  a..  E.  cos  [(1  -  n.)t  +  0-1 

2  ij  j  r  y  ]J 

=  a..  E.  f cos  f(l  -  fi.  -  fi. )  t  +  0.  +  0. 1  cos (£2.  t  -  0. ) 

2  ij  jl  T  ij  i  j]  '  i  i 

-  sin  |(1  -  £2i  -  £2^  )t  +  0.  +  0  J  sin(fi.  t  -  6^)  J 

(A -5) 


and  a  similar  equation  with  the  subscripts  i ,  j  interchanged.  Then  by  equating  the 
coefficients  of  cos  ( ftf  t  -  0  )  ,  sin  ( ft  t  -  0  )  ,  r  =  i ,  j  we  get 


2 ft.  0 . 

E.  = 

a 

r\ 

a.. 

E. 

cos 

(1 

-  ft.  - 

ft.)t 

+ 

0. 
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0. 
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2ft. 
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ij 
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2ft .  0 . 
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ft.  -  ft 
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0. 
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J  J 
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+ 
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l 
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JJ 

(A-6) 


♦Struble,  R.A.,  Nonlinear  Differential  Equations,  New  York, McGraw-Hill  Book 
Company,  1962,  Chapters,  pp.  221-227 
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It  can  be  easily  shown  that 


ft.  a..  E?  =  ft.  a..  E2 
i  Ji  i  J  U)  J 


e.  =  0. 

1  J 


Substituting  into  Eq.  (A -6),  it  becomes 

1/2 


2Vi  =  f(-iLnp)  cos[(l  -  -  vt+  2ei\ 


a. .  a.,  ft. 

“  1J  ‘ 

> 

“j 


1/2 


(a. .  .  a. . .  4  .  \ 

e.  cos  (1  -  ft.  -  ft.)t  +  20. 
ft.  i  r  i  j7  ij 


(A— 7 ) 


(A-8) 


The  values  of  ft.  ,  ft.  ,  a. . ,  a..  ,  a  will  determine  whether  E.,E.  in  the  above  equation, 
i  J  ij  Ji  i  J  M 

has  unbounded  solutions  or  not.  We  note  that  Eq.  (A-8)  has  exactly  the  same  form  as 

that  obtained  bv  Struble  for  a  single  Mathieu’s  equation.  If  a.,  a.,  ft. /ft.  is  less  than  0 

ij  Ji  i  J 

from  Struble1  s  result,  all  solutions  of  E.  are  bounded;  if  a.. a.,  ft.  ft.  is  greater  than  0, 

i  ij  ji  i  j 

its  stability  criterion  is 


1  -  ft.  -  ft.| 

i 


a  | 

/ a.,  a..  \ 

1  U...J 1 

s 

2  ' 

V 

Stable 

Unstable 


(A-9) 
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